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ABSTRACT 

PLh \ We present new lightcurves of the massive hot Jupiter system WASP-18 obtained 

r"l"1 ■ with the Spitzer spacecraft covering the entire orbit at 3.6 /jm and 4.5 yum. These 

lightcurves are used to measure the amplitude, shape and phase of the thermal phase 
effect for WASP-18 b. We find that our results for the thermal phase effect are limited 
to an accuracy of about 0.01 per cent by systematic noise sources of unknown origin. 
At this level of accuracy we find that the thermal phase effect has a peak-to-peak 
amplitude approximately equal to the secondary eclipse depth, has a sinusoidal shape 
and that the maximum brightness occurs at the same phase as mid-occultation to 
within about 5 degrees at 3.6 /jm and to within about 10 degrees at 4.5 /xm. The shape 
and amplitude of the thermal phase curve imply very low levels of heat redistribution 
within the atmosphere of the planet. We also perform a separate analysis to determine 
the system geometry by fitting a lightcurve model to the data covering the occultation 
and the transit. The secondary eclipse depths we measure at 3.6 /im and 4.5 /im are in 
good agreement with previous measurements and imply a very low albedo for WASP- 
18 b. The parameters of the system (masses, radii, etc.) derived from our analysis are 
in also good agreement with those from previous studies, but with improved precision. 
We use new high-resolution imaging and published limits on the rate of change of the 
mean radial velocity to check for the presence of any faint companion stars that may 
affect our results. We find that there is unlikely to be any significant contribution 
to the flux at Spitzer wavelengths from a stellar companion to WASP-18. We find 
that there is no evidence for variations in the times of eclipse from a linear ephemeris 
greater than about 100 seconds over 3 years. 
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1 INTRODUCTION 

Hot Jupiters are currently at the forefront of observa- 
tional studies that can provide meaningful tests for mod- 
els of exoplanet atmospheres. The atmospheric tempera- 
tures for a typical hot Jupiter orbiting a solar-type star 
with a period of 3 days can be up to 1500 K. For tran- 
siting hot Jupiters this makes it feasible to measure the 
planet-star flux ratio directly from the depth of secondary 
eclipse in the lightcurve due to the occultation of the exo- 
planet by the host star. Early results with the Spitzer Space 
Telescope confirmed the existen ce of secondary ec lipses 
in the li ghtcurves of HP 209458 (jDeming et al.ll2005T ) and 
TrES-1 (|Charbonneau et alj|2005h with the expected depth 
~ 0.5 per cent at mid-infrared wavelengths. The secondary 
eclipse depth has now been measured using Spitzer for more 



than 20 hot Jupiters jCowan fc Agol|[2uTll ). Comparison of 
these observations with atmospheric models has been used 
to reveal the diversity of h ot Jupiter atmospheres wit h re- 
gard to their composition (|Madhusudhan et al.ll201lT ). the 
presence or absence of a tempe rature inversion in their atmo- 
spheres, (|Knutson et aljfeoiol). and their albe dos and heat 
recirculation efficiencies 1 Cowan fc Ago]|[201ll ). 

The secondary eclipse depth for a hot Jupiter at in- 
frared wavelengths measures the brightness temperature of 
the hemisphere facing the star - the "day-side" - integrated 
over the visible hemisphere. This brightness temperature, 
Tday, will depend on the pattern of emission over the day 
side, the spectral energy distribution (SED) of this emission, 
the Bond Albedo, At, and the efficiency with which heat is 
redistributed to the night-side of the planet. Observations at 
several wavelengths, particularly near-infrared observations 
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near the peak of the day-side SED, reduce the extent to 
which we must rely on models to account for the conversion 
from brightness temperature to effective temperature whe n 
interpreting these observations l|Madhusudhan et al.ll20l"ll ). 
In general, it is not possible to disentangle the degeneracy 
between Bond albedo and heat recirculation efficiencies from 
secondary eclipse observations alone. One exception is the 
case of HD 189733, in which very high quality Spitzer ob- 
servations of the secondary eclipse reveal asymmetries that 
can be inverted to produce a map of th e brightness temper- 
ature on the day s ide of HD 189733 b |Maieau et alj|2012l ; 
Ide Wit et alj|20l3 ). Apart from this exceptional case, it is 
currently only possible to obtain information on the redistri- 
bution of heat in the planet's atmosphere by observing the 
thermal phase effect - the variation in infrared brightness of 
the system as a function of orbital phase. 

There are several different ways to parametrize the re- 
distrib ution of heat from the da y-side to the night-side of a 
planet dSpiegel fc Burrowsil2010h . Here we use the parameter 
P n , the fraction of the incident energy that is transported 
to the night side of the planet. Plausible values of this pa- 
rameter vary from up to 0.5. A value of P n — would 
imply a night-side brightness temperature T n i g ht <C Tday, 
whereas P n = 0.5 implies T lligh t ~ TdayQ If Pn > then this 
suggests that winds at some level in the atmosphere move 
heat around the planet. In practise, very high efficiencies for 
strongly irradiated planets are unlikely because the winds 
that transport heat to the night-side will dissipate some o f 
their energy through turbulence or shocks |Goodmanl l2009) ■ 
Nevertheless, some redistribution of energy from the day- 
side to the night-side is likely, and may lead to significant 
offsets between the sub-stellar point and the h ottest regions 
of the atmosphere l|Cooper fe Showmanll2005T ). This will be 
observed in the thermal phase effect as an offset in the phase 
of maximum brig h tness from opposition. 

ICowan et alj l|2007t) used 8 separate observations with 
Spitzer spread throughout the orbit of the three hot Jupiter 
systems HD 209458, HD 179949 and 51 Peg to measure their 
thermal phase effect. They were able to place useful upper 
limits on the phase variation in 51 Peg and HD 209458 and 
to detect a variation with a peak-to-trough amplitude of 
0.14 per cent in HD 179949. HD 179949 is a non-transiting 
hot Jupiter, so the inclination of the orbit and the radius of 
the planet are unknown, but even allowing for this uncer- 
tainty, the observed amplitude of the phase variation pro- 
vided an upper limit of P n < 0.21 and shows that the hottest 
point is near the sub-stellar point. HD 209458 is a transit- 
ing system, so the inclination of the orbit and the rad ius 
of the exoplanet are known. This allowed ICowan et al.l to 
translate the upper limit on the amplitude of the thermal 
phase into a lower limit P n > 0.32, and thus establish that 
apparently different hot Jupiters are likely to have a variety 
of Pn values. 

lHarrington et all (|2006t ) detected the phase variation of 
the planet v Andb using the MIPS instrument on Spitzer at 
24/im. With additional data and an improved understand- 
ing of the systematic noise sources in MIPS, they were able 
to refine their estimate of the amplitude of the phase vari- 



The brightness temperatures are not necessarily equal since the 
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ation and show that there is a large (~ 80°) phase off- 
set bet ween the time of ma ximum brightness and oppo- 
sition (|Crossfield et al.l [2010). The inclination and radius 
of v Andb are not known accurately because it is a non- 
transiting exoplanet. This can result in a large uncertainty 
in the value of P n inferred from the amplitude of the phase 
curve |Burrows et alj [20081 ). However, the large phase off- 
set observed in v And b implies a large P n despite the large 
amplitude of the phase variation. 

The thermal phase curve of HD189733 b has been ob- 
serve d using Spitz er at 3.6/im an d 4.5/i m l|Knutson et al.l 
l2012h. at 8um jKnutson et all 120071 ) and at 24/xm 
jKnutson et~afl l2009h . There are also multiple observa- 
tions of the tran sits and eclipses at 8/j.m for this system 
dAgol et alj|20ld) . The combined analysis of these results 
bv lKnutson et all l|2012t ) shows that heat recirculation from 
the day-side to the night-side is efficient for this relatively 
cool hot Jupiter (Td ay ~ 1200 K) and that this recircula- 
tion leads to a peak in the thermal phase effect that occurs 
~ 25 degrees before opposition. 

ICowan et all (|2012h obtained Warm Spitzer photom- 
etry covering the complete 26 hour orbit of the very hot 
Jupiter WASP-12 b at 3.6 /im and 4.5 /im. They found that 
their interpretation of the lightcurves depends on the as- 
sumptions made about the nature of the systematic noise in 
the lightcurves and that red noise is the dominant source 
of uncertainty in their analysis. Nevertheless, they were 
able to show that the thermal phase variation in WASP- 
12 b is large, indicative of poor day-to-night heat redistri- 
bution (P n < 0.1). The small offset they observe between 
the phase of maximum brightness and secondary eclipse 
(16 ± 4 degrees) in the 4.5 ^m data is consistent with this 
interpretation. The phase offset at 3.6 /im could not be de- 
termined unambiguously from their data. 

Although thermal phase curves are only currently avail- 
able for a few systems, there does appear to be a pattern of 
weak recirculation for the hottest planets. I Cowan fc Agoll 
l|2011f ) have looked for trends in the value of Td ay /T su b, 
where T au b is the equilibrium temperature of the sub-stellar 
point, in a sample of 24 transiting exoplanets with secondary 
eclipse measurements. This quantity will depend on both the 
albedo of the planet and the recirculation efficiency, but a 
large value can only be obtai n if both the alb edo and the re- 
circulation efficiency are low. ICowan fc Agoll found that this 
is the case for all of the 6 hottest planets (Td ay ^ 2400 K) in 
their sample and point out that this is, in general terms, 
the expected behviour given that the radiative timescale 
scales as T~ 3 whereas the advective timescale (which they 
assume to be of-order the local sound speed) scales as T -0 ' 5 . 
This simple scaling argument does not explain the apparent 
transi tion in behaviour at Td ay « 2400 K but iPerna et al.l 
(2012) do observe a transition at about this temperature in 
their suite of three-dimensional circulation models for hot 
Jupiters. This is mainly due to the change in the ratio of 
the radiative and advective timescales, with the presence of 
an atmospheric inversion playing a lesser role in determining 
the recirculation efficiency. This transition may also be re- 
lated to the onset of ionisation of alkali metals in the planet's 
atmo sphere, leading to severe magnetic drag (|Perna et al.l 
2010). There is an ongoing debate as to whether the result- 
ing Ohmic dissipation can transport energy into the interior 
of the planet and so explain the very large radii observed for 
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some hot Jupiters ([Laughlin et al.ll201ll;lRauscher fc Menoul 
l2012l ; |Perna et al J 12012 : iHuang fe Cumminell2012l ). 

Here we present Warm Spitzer photometry covering the 
complete orbit of the very hot Jupiter WASP-18 b. This ex- 
oplanet is unusual for its combination of short orbital pe- 
riod (0.945 d) and high mass (10M,j up ), which results in 
strong tidal interac ti ons between the planet and the star 
|Hellier et all 120091 ). ISouthworth et all l|2009h derived ac- 
curate masses and radii for the star and planet in the 
WASP-18 system based on high quality optic al photometry 
of the transit and the spectroscopic orbit from lHellier et all . 
iNvmever et~aH i|201ll ) used Spitzer photometry in all four 
IRAC bands covering the secondary eclipse of WASP-18 to 
measure the brightness temperature of the day side from 
3.6 fim to 8.0 fim. The high brightness temperatures derived 
(~ 3200 K) imply that WASP-18 b has near-zero albedo and 
almost no redistribution of energy from the day side to the 
night side of the planet. 

Our primary aim is to use our Warm Spitzer photome- 
try at 3.6 /im and 4.5 /im to measure the amplitude, phase 
and shape of the thermal phase effect. We also use these 
data to re-measure the secondary eclipse d epths at 3.6 fim 
and 4. 5 fim for comparison with the results of lNvmever et all 
(2011). We measure the times of eclipse and transit from our 
data and from new optical photometry of several transits 
and use these together with published times of mid-eclipse 
to re-measure the eclipse ephemeris and to look for possi- 
ble variations in the period. We consider the likelihood that 
WASP-18 has a stellar companion based on published radial 
velocity data and new high-resolution imaging at H-band 
and K-band. The contamination of the lightcurve for a hot 
Jupiter system by a companion star has the potential to bias 
the results obtained if not properly accounted for. Compan- 
ions stars may also play a role in t he formation and evo- 
lution of some hot J upiter systems (|Fabrvckv fe Tremaind 
l2007l ; lMardiiii3l2007l ). Our analysis also provides an accu- 
rate characterisation of the primary eclipse (transit) which 
can be used in combination with other data to re-measure 
the mass and radius of the star and planet. 



2 OBSERVATIONS 
2.1 Spitzer photometry 

We were awarded Spitzer General Observer time during Cy- 
cle t Qto observe tw o complete orbits of WASP-18 with IRAC 
l|Fazio et al.l Il998l ). one orbit with each of the two IRAC 
channels operating during the warm mission. Observations 
with channel 1 (3.6/im) were obtained on 2010 January 23, 
observations with channel 2 (4.5/im) were obtained on 2010 
August 23. On both dates, 243 200 images with an expo- 
sure time of 0.36s were obtained in sub-array mode. The 
total duration of each sequence of observations is 29 hours. 
In addition, a sequence of 64 sub-array images also with an 
exposure time of 0.36s were obtained immediately after the 
observations of WASP-18 at a slightly offset position. These 
were used to check for hot pixels or other image artifacts on 



the detector. We used Basic Calibrated Data (BCD) pro- 
cessed with version S18.18 of the Spitzer IRAC pipeline for 
our analysis. 

2.2 AO imaging 

We obtained adaptive optics high resolution H- and K-hanA 
images of WASP-18 using the NICI instrument at Gemini 
South. The observations were obtained on the night of 2010 
December 27 under good seeing (0.5 -0.6 arcsec). The in- 
strument was configured with the CLEAR focal plane mask 
and the H 50/50 beam-splitter, and we used the narrow band 
filters iCcont (2.2718 fim) and Fell (1.644 fim) in the red and 
blue channels, respectively. We observed the target at five 
dither positions corresponding to the center and corners of a 
square of side 6 arcsec. At each dither position we obtained 
three images consisting of the co-addition of 3 exposures 
of 1.5 s. The full-width at half-maximum (FWHM) of the 
point spread function (PSF) in these images was 0.065 arcsec 
in the Fell filter and 0.073 arcsec in the K con t filter. 

Data reduction consisted of subtracting a sky image, di- 
viding by the flat field, and fixing bad pixels by interpolating 
from neighbouring pixels. The sky image was created from 
the median combination of the images at the different dither 
positions after masking out the star signal in each image. 
The reduced images were registered to a common position 
and field orientation and then combined using the median 
value of each pixel. No other point source was detected in 
the resulting images. The sensitivity of our AO imaging to 
detect faint companions was determined by first computing 
the median absolute deviation of the pixel values within an- 
nuli of various radii and width equal to one PSF FWHM. 
The resulting contrast curve was then properly scaled, and 
verified to be adequate, by adding and recovering (by vi- 
sual inspection) fake companions in the images at various 
separation and with various contrasts. In doing this last ex- 
ercise we used both the K- and //-band images to differenti- 
ate speckles from true companions, which display a different 
chromatic behaviour. Using this approach, we estimate the 
detection limits in difference of magnitudes to be 4.0 mag at 
^0.2", 5.4 mag at ^0.4" and 6.0 mag at ^0.5". Finally, we 
note the presence of a faint ghost in the image at 13 pixel 
(0.23") separation and contrast of ~4.1 mag in the iCcont 
filter and ~5.9 mag in the Fell filter. 

2.3 TRAPPIST photometry 

Five transits of WASP-18 b were o bserved with the 
60 cm robotic tel escope TRAPPISlEI ( Jehin et al. 

I 1201 ll : 

iGillon et all |2012| ) located at ESO La Silla Observatory 
(Chile) . TRAPPIST is equipped with a thermo-electrically- 
cooled 2k x 2k CCD camera with a field of view of 22' x 22' 
(pixel scale = 0.65 arcsec pixel -1 ). The first four transits 
were observed in an Astrodon 'J + z' filter that has a trans- 
mittance > 90 per cent from 750 nm to beyond 1100 nm, the 
red end of the effective bandpass being denned by the spec- 
tral response of the CCD. The last transit was observed in 
the Sloan z' filter. For all transits, the telescope was slightly 
defocused to minimize pixel-to-pixel effects and to optimize 
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Table 1. Details of the transit lightcurves obtained with TRAP- 
PIST for WASP-18. For each lightcurve we list the observation 
date, the filter used, the number of measurements, and the expo- 
sure time. 



Date 


Filter 


N p 




2010 Sep 30 


I+z 


712 


12 


2010 Oct 02 


I+z 


977 


8 


2010 Dec 23 


I+z 


688 


6 


2011 Jan 08 


I+z 


648 


6 


2011 Nov 10 


z' 


624 


10 



the observational efficiency. For each run, the stellar images 
were kept on the same pixels, thanks to 'software guiding' 
system deriving regularly astrometric solutions on the sci- 
ence images and sending pointing corrections to the mount 
if needed. Table [T] gives the logs of these TRAPPIST obser- 
vations. 

After a standard pre-reduction (bias, dark, flatfield cor- 
rection), the stellar fluxes were extracted from the images 
using the ira f/daophot3 aperture photometry software 
i|Stetsod[l"987l ). For each transit, several sets of reduction 
parameters were tested, and we kept the one giving the 
most precise photometry for the stars of similar brightness 
as WASP-18. After a careful selection of reference stars, dif- 
ferential photometry was obtained. 

3 SPITZER DATA ANALYSIS 

3.1 Conversion to flux and noise model 

We converted the BCD images from units of MJy/steradian 
to mjy using the values for the pixel size at the centre of 
the subarray provided in the image headers (1.225" x 1.236" 
for channel 1, 1.205" x 1.228" for channel 2). We used the 
raw images together with the values of the gain and readout 
noise for each channel to calculate the noise level in each 
pixel assuming Poisson counting statistics. 

3.2 Image times 

BCD data in sub-array mode are delivered as FITS files con- 
taining a data cube of 64 images of 32 x 32 pixels per file. 
We used the FITS header keyword BMJD.OBS to assign a 
Barycentric UTC modified Julian date (BMJD) to the start 
of the exposure for the first image in the data cube. The 
BMJD of the mid-exposure time for each image in the data 
cube was then calculated using the values for the start and 
end times of the integration from the FITS header (AINT- 
BEG and ATIMEEND) to calculate the time taken to ob- 
tain the 64 images and assuming that these images were 
uniformly spaced in time. 

Long observing sequences such as the ones we have used 

4 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



for WASP-18 cannot be executed using standard observing 
modes so multiple instrument engineering requests (IERs) 
are used to obtain the data. The images obtained in the sec- 
ond of the two IERs used for our WASP-18 observations do 
not have the coordinates of the target in the FITS header 
and so the light-time correction to the solar system barycen- 
tre are incorrect for these images. We calculated the light- 
time correction for the images obtained before the interrup- 
tion from the difference in the keyword values BMJD_OBS — 
MJD_OBS. We then use a linear extrapolation of this light- 
time correction as a function of MJD_OBS to calculate the 
BMJD of the images obtained after the interruption based 
on their MJD_OBS values. The uncertainty in the exposure 
time introduced by this procedure is negligible. 



3.3 Outlier rejection 

We compared each image to the other 63 images in the same 
data cube in order to identify discrepant data points in the 
images. We are particularly concerned here with identifying 
discrepant pixel values that may affect the photometry of 
the target. As the target moves during the sequence of 64 
images, we use a robust linear fit (least absolute deviation) 
to the 64 pixel values from each file for each pixel to predict 
the expected value for each pixel value in each image. We 
then flag the pixels in each image that deviate from their 
expected value by more than 5 times their standard error. 
We find that the number of pixels flagged using this method 
is much larger than expected given the known incidence of 
cosmic ray hits on the IRAC detectors. This discrepancy 
is due to a few pixels well away from the target position 
that are noisier than predicted by our noise model. As these 
pixels have no effect on our photometry and a negligible 
effect on the estimate of the background level we ignore this 
discrepancy. We also flagged any pixels in our images that 
are flagged as bad pixels in the "Imask" file provided for 
each BCD file by the Spitzer IRAC pipeline. 



3.4 Sky background estimate 

We use the mean of the image pixel values excluding those 
within 10 pixels of the target position to estimate the back- 
ground value in each image. Values more than 4 standard 
deviations from the mean and flagged pixels were ignored 
in the calculation. We used a Gaussian fit to a histogram of 
these pixel values to estimate the standard deviation of the 
background pixel values, (Jb g . The number of points used to 
estimate the background was ~ 700. Typical values of o~b g 
are 0.0033 mjy/pixel for channel 1 and 0.0025 mjy/pixel for 
channel 2. 



3.5 Aperture photometry 

We tried three different methods to measure the location of 
the star on the detector, the daophot cntrd and gcntrd 
algorithms and a least-squares fit of a bivariate Gaussian 
distribution to an 11x11 sub-image centred on the nomi- 
nal star position. We refer to this latter algorithm as the 
GAUSS2D method. We used a fixed value for the full-width 
at half-maximum of FWHM= 1.25 pixels for both axes of 
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Channel 1 before 



Channel 2, before 




Figure 1. Images obtained before and after our WASP-18 obser- 
vations. All images are linearly scaled (inverse grey scale) between 
and 10 MJy/sr. The "before" images are 30 s exposures in the 
region of IC 2560. The "after" images are the median of the 64 
offset sky images with an exposure time of 30 s each. 



the Gaussian profile in the GAUSS2D method based on the re- 
sults of fitting the images with the FWHM as a free param- 
eter. The CNTRD algorithm determines the position where 
the derivatives of the image values go to zero. The GCNTRD 
algorithm fits a Gaussian profile to the marginal x and y dis- 
tributions of the image values. We set the input parameters 
to CNTRD and GCNTRD to run on a sub-image of 5 x 5 pixels 
around the target position. We compare the performance of 
these different algorithms below. 

We used the IDL Astronomy Users li brar^fl implem en- 
tation of the daophot aper procedure (|Stetsonlll987l ) to 
perform synthetic aperture photometry on our images. We 
used the 2006 November version of this procedure which al- 
lowed us to use the option to use an exact calculation of the 
intersection between a circular aperture and square pixel 
for correct weighting of pixels at the edge of the aperture. 
We used 13 aperture radii uniformly distributed from 1.5 
pixels to 4.5 pixels. The results we obtained for larger aper- 
ture radii were not useful because the lightcurves have much 
lower signal-to-noise due to the additional background noise. 
Fluxes measured from images containing any flagged pixels 
in the aperture were rejected from further analysis, although 
much less than 1 per cent of images were affected in this way. 
The median number of rejected pixels per image is 2. 



3.6 Image persistence 

Our IRAC images are affected by image persistence, partic- 
ularly the channel 1 images. This can be seen in the offset 
sky images obtained immediately after our WASP-18 ob- 
servations (Fig. [TJ. The resulting image artifacts are more 
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diffuse than the stellar images and look more like the loga- 
rithm of the point spread function, as described in the IRAC 
instrument handbook (Version 2.0.1, p. 116) . The image ar- 
tifact in the channel 1 offset sky images has up to about 0.6 
mjy per pixel in the channel 1 image and a total of about 
6 mjy within an aperture with a radius of 5 pixels. For com- 
parison, a typical channel 1 image of WASP-18 has a peak 
flux of 50-70mJy/pixel and a total flux of 163 ± 3mJy 
within an aperture of the same radius, so the image artifact 
affects the photometry of WASP-18 by a few percent. For 
the channel 2 data the corresponding figures are a total of 
about 1 mjy/pixel in the artifact compared to 42mJy/pixel 
in the peak and a total flux of 103 ± 1 mjy in the images of 
WASP-18 so for this channel image persistence affects the 
photometry by about 1 per cent. 

The IRAC instrument handbook describes the be- 
haviour of image persistence artifacts in the IRAC arrays 
during the warm mission. In channel 1 the artifacts decay 
exponentially with a timescale of about 4.5 hours. Channel 
2 residuals start out as positive, but then become negative 
with a decay timescale of a few minutes. 

For channel 1 data we make the assumption that the 
image artifact will be approximately constant after some 
time during the exposure sequence comparable to the de- 
cay timescale. To correct for the effect of image persistence 
we create a "master offset image" from the median of the 
64 offset sky images, subtract the background value from 
this image and then subtract the result from the images of 
WASP-18. This correction will be inaccurate for some frac- 
tion of the data at the start of the observing sequence while 
the image persistence builds-up. We discuss this point fur- 
ther below. 

For the channel 2 data it is not clear how the image ar- 
tifact affects the photometry of WASP-18. The interval be- 
tween the end of the observing sequence for WASP-18 and 
the start of the offset sky image is 49 s, which is comparable 
to the decay timescale for the artifact. The exact form and 
timescale for the decay of the artifact is not known so it is not 
even possible to make a precise estimate of contribution of 
the image artifact to the measured flux for WASP-18. How- 
ever, it is likely that the image artifact contributes less than 
2 per cent given the decay timescale for this feature is a few 
minutes. We did attempt to measure the decay timescale 
from the data taken subsequent to our own observations, 
but the artifact was not detectable in those images. For the 
channel 2 data we treat the contribution of the image persis- 
tence artifact as an additional source of uncertainty in our 
analysis. 

3.7 Initial assessment of the data 

The flux of WASP-18 measured with an aperture radius of 
4.5 pixels is shown for both channels in Fig. [2] Also shown 
in this figure are the positions of the star on the array cal- 
culated using the GAUSS2D method. The coordinates x and 
y are measured relative to the centre of a corner pixel in 
the sub-array. The form of the variations in the x,y posi- 
tions measured using the GCNTRD and CNTRD algorithms 
are similar, but the amplitudes of the variations are less and 
there is an offset between these values and the results of the 
GAUSS2D method. For example, the y positions measured 
for the channel 1 data using the GCNTRD method have me- 
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dian value of 14.95 with 98 per cent of the data in the range 
y — 14.78-14.95, cf. a median value of 14.88 and range 
y = 14.65-14.99 for the GAUSS2D method. 

The feature that stands out from Fig. [2] is the well- 
known correlation between the measured flux and the posi- 
tion of the star on the detector, particularly in channel 1. 
This position-dependent sensitivity variation (PDSV) makes 
it difficult to see the transit and secondary eclipse in these 
"raw" aperture flux measurements. The channel 2 data ap- 
pear to be less affected by PDSV, so the transit and sec- 
ondary eclipse can be seen in the raw flux measurements. 
PDSV is a combination of the "pixel phase effect" described 
in the IRAC instrument handbook (Version 2.0.1, p. 45) and 
pixelation noise. The pixel phase effect is a variation in the 
sensitivity of each detector pixel that depends on the dis- 
tance of the stellar image from the centre of the pixel. The 
motion in the x and y directions for our channel 2 data result 
in a smaller variation in the distance of the star from the cen- 
tre of the pixel compared to the channel 1 data, which may 
part ly explain why the d ata quality is better in this chan- 
nel (|Anderson et al.ll201lf ). Pixelation noise affects synthetic 
aperture photometry with small aperture radii because for 
a pixel at the edge of the aperture, the fraction of the flux 
falling within that pixel is not the same as the fraction of 
the aperture within the pixel. 

3.8 Analysis of the thermal phase effect 

3. 8. 1 Correction for PDSV 

The usual method developed to correct for PDSV in IRAC 
data for observations of the secondary eclipses of hot 
Jupiters is to include parameters in the least-squares fit of an 
eclipse model to the data to represent the PDSV. This is usu- 
ally a simple linear or quadratic relation between sensitivity 
and each of the coord i nates x and y ( e .g.. | Beerer et al.l[201l] : 
lAnderson et alj|201ll ). iBallard et ail (|201(t ) have developed 
an alternative method to correct for PDSV in their Warm 
Spitzer 4.5/im observations of GJ436. They created a pixel 
sensitivity map from the data themselves. This approach was 
straightforward in the case of GJ 436 because the lightcurve 
of the target is expected to be constant apart from the possi- 
ble presence of transits affecting a smal l fraction of th e data. 
The pixel sensitivity map generated bv lBallard et al.l for the 
IRAC channel 2 shows complex structure that they describe 
as "corrugation . . . low-level sinusoidal-like variations with 
a separation of approximately 5/100ths of a pixel between 
peaks". A similar concept based on bi-linear interpolation 
rat her than a smoothed look-up table has been developed 
bv IStevenson et "all ((2012h and applied to Spitzer photome- 
try of HD 149026 b. 

For our WASP- 18 data we are interested in character- 
ising the amplitude and shape of the phase variation as well 
as measuring the shape and depths of the transit and sec- 
ondary eclipse. The phase variation has a period comparable 
to the length of the observing sequence so it is important 
to understand any correlations between the correction for 
the PDSV and the parameters of the lightcurve model. The 
phase variation of WASP-18 can be modelled approximately 
as a sinusoidal variation in flux with P=0.94d. For exam- 
ple, in the worst case scenario, if the position of the target 
on the detector also varied approximately sinusoidally with 



approximately the same period, then it would become im- 
possible to determine whether any variation in the measured 
flux with P « 0.94 is due to the flux variation of WASP-18 
or due to the PDSV. 

Understanding the correlations between the lightcurve 
model parameters and the correction for PDSV is prob- 
lematic in the case of our WASP-18 data because of the 
large number of parameters required to model the complete 
lightcurve plus the large number of parameters that may 
be required to characterise the complex structure ("corru- 
gation") in the PDSV. One method we expe rimented with 
was to use the pixel sensitivity map method of lBallard et alj 
(2010) applied to the residuals to a least-squares fit of a 
lightcurve model. This approach can be applied iteratively 
until the solution and pixel sensitivity map converge. The 
problem with this approach is that it becomes difficult to 
identify correlations between the lightcurve model parame- 
ters and the pixel sensitivity map. We avoided this problem 
by excluding the data during the transit and occultation 
from our analysis of the thermal phase effect. The main ad- 
vantage of this approach is that fitting a model to the data 
between the eclipses can be reduced to a linear least-squares 
problem. This makes is straightforward to find the best so- 
lution of the problem and to investigate the correlations be- 
tween the free parameters of the model. 

3.8.2 Model for PDSV and the thermal phase effect 

Our model for the measured magnitude of the system be- 
tween the eclipse and transit is 

mi = c ,o + Ylj=i a J c°s{j<f>i) + YlkS bk sm (fc0O 

+ E£i «K) + ESr co, K p K (yl) (1) 

where rrij is the magnitude of WASP-18 at time ti\ <f>i = 
2n(ti — To)/P is the orbital phase relative to the time of 
mid-transit, To; p„ is a Legendre polynomial of order n; 
x'i = (xi - x)/(x max - :r m in) and similarly for y[ (x min is 
the minimum value of x it etc.). We use the values To = 
BMJD 54220.98163 for the time of m id-transit and P = 
0.94145299 d for the orbital period from lHellier et all (|2009T) 
throughout this paper unless otherwise stated. By using Leg- 
endre polynomials and normalized coordinates (x'i,y'i) we 
find that we can use singular value decomposition to find 
solutions of this least-squares problem for Legendre polyno- 
mials up to at least 12 th order. This is sufficien t to model the 
corrugations with a scale of 0.05 pixels seen bv lBallard et al.l 
if they are present in our data. 

3.8.3 Linear decorrelation against position 

In Fig. Owe show the result of using the simplest reasonable 
model for our data, in which the magnitude of the phase 
variation varies sinusoidally and the PDSV is linear in x 
and y, i.e., 

m t = c ,o + ai cos(0,) + bi sin(0 t ) + c lfi p 1 (x[) + c ,ipi (?/•)• 

Note that our calculations are done using magnitudes, but 
we plot the results as fluxes and quote parameter values 
in unit of percent. The least-squares fit of this model to 
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Figure 2. The flux of WASP-18 measured in IRAC channel 1 (left panel) and channel 2 (right panel) measured with a circular aperture 
of radius of 4.5 pixels.. The position of the star on the array measured using the GAUSS2D method is shown below each panel. For clarity, 
we have only plotted a random selection of 1 per cent of the data here. Dashed lines indicate the start and end times of the transit and 
secondary eclipse assuming a duration of 0.08 d for each and assuming that the secondary eclipse occurs at phase 0.5. 



the unbinned aperture photometry outside of eclipse and 
transit is used to determine the coefficients of the model for 
the PDSV. We then apply this correction to all the data. 
The results in Fig. [3] are for an aperture radius of 2.5 pixels 
for both apertures and target positions measured using the 
GAUSS2D method. This is the combination of aperture radius 
and positions that gave the lowest RMS residual for the data 
between the transit and eclipse. Results for other apertures 
and for CNTRD and GCNTRD methods are similar. It is clear 
that a linear correction is insufficient to fully remove the 
effect of the PDSV, but this simple model does show clearly 
some features of our data. Firstly, we see that the eclipse and 
transit are clearly visible in both channels. Two transits are 
visible in the channel 2 data but the first transit is not seen 
clearly in the channel 1 data because there is a large "ramp" 
affecting the first few hours of the data. There is a cosine- 
like variation in flux observed in both channels with the 
maximum flux occurring near phase 0.5 (secondary eclipse). 
Part of this signal is the phase variation we wish to measure. 
However, there must also be some instrumental component 
or other systematic effect that contributes to this variation 
because the phase variation due to the planet cannot have 
an amplitude larger than the secondary eclipse depth. 



The obvious suspect for the systematic noise source in 
the channel 1 data is the image artifact shown in Fig. [T] 
The "ramp" is the right size (w 3 per cent) and builds up 
over the same sort of timescale as the known decay timescale 
of this artifact. Our interpretation of this lightcurve is that 
the image artifact builds up over the first 6-8 hours before 
reaching an approximate equilibrium between the arrival of 
new photons from WASP-18 and its own decay timescale. 

We have tried several methods to account for this ramp- 
like feature in the data but none of these methods is any 
better than the more pragmatic approach of simply exclud- 
ing the first 6-8 hours of the channel 1 data. Without go- 
ing into the details of these various methods, we can state 
here that we almost always found that the amplitude of the 
phase variation measures in channel 1 was similar to the 
depth of the secondary eclipse and often was slightly larger. 
It is possible to create models for the systematic noise in 
the channel 1 lightcurves that achieve more physically real- 
istic (lower) values for the amplitude of the phase variation, 
but these models are not based on any physical model of 
the instrumental noise, i.e., they are arbitrary, and they re- 
quire several additional free parameters that are often not 
well constrained by the data or any physical understand- 
ing of what these parameters represent. The overall quality 
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Figure 3. The flux of WASP-18 measured in IRAC channel 1 (left panel) and channel 2 (right panel) after a linear correction for position 
dependent sensitivity variations and a simple sinusoidal model for the phase variation. The zero-point of the flux scale is set from the 
mean flux during secondary eclipse. Observations obtained during the transit and secondary eclipse (small points) were excluded from 
the calculation of the coefficients for the decorrelation. The data have been binned into 0.0025 d bins for display purposes only. 



of the decorrelated lightcurve obtained with these arbitrary 
and complex models of the the instrumental noise is also not 
much better than the best results presented below for the 
partial lightcurve. Clearly, a more complete understanding 
of the instrumental noise in IRAC for warm mission obser- 
vations would be a great help for the interpretation of our 
data, but in the absence of this we present the results for 
the partial lightcurve and make an attempt to quantify the 
extent to which instrumental noise introduces systematic er- 
rors in our results. 



3.8.4 Optimum decorrelation against position. 

We used the model given in equation (1) to fit the data 
for channel 1 and channel 2 excluding data within 0.05 d 
of mid-transit and mid-eclipse and also excluding the first 
60 000 observations (7.3 hours of data) for channel 1. We 
used N cos = 2 and N B in = 1 to model the phase variation 
of WASP-18. The sine term allows for a phase shift from 
phase 0.5 for the time of maximum brightness and the first 
harmonic of the cosine variations (coefficient a-z) allows for 
some optimisation of the shape of this phase variation. To 
model the PDSV we tried N x = N xy = AT y /2 = 1, 2, ... 6. 
We use N y = 2N X because there is a larger range of motion 



in the y direction. We fit these models to the lightcurves for 
all combinations of aperture size and position measurement 
methods. We used the Bayesian information criterion (BIC) 
to identify the combination of (N x , N y , N xy ) that provides 
the best compromise between quality of fit and number of 
free parameters for a given lightcurve. We calculated the 
BIC using the expression 

BIC = x 2 + A r P a r log e (A), 

where N pal is the number of free parameters and N is the 
number of observations. We used the RMS of the residuals to 
identify the aperture size and position measurement method 
that give the best lightcurves. For both channels we find the 
best results are obtained for (N x , N y , N xy ) — (5, 10, 5) with 
positions measured using the CNTRD method and an aper- 
ture radius of 2.25 pixels. These models and lightcurves are 
shown in Fig. [4] and the parameters of interest are given in 
Table The standard error estimates given in Table [2] ac- 
count for the correlations between parameters (|Press et al.l 
Il992l ) although the correlation coefficients between the pa- 
rameters in this table and all other parameters in the model 
are small (< 0.3). 

In Fig. [5] we show how the parameters ai , 02 and 61 
obtained for (N x , N y ,N xy ) — (5, 10, 5) vary as a function of 
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Table 2. Results of linear least squares fits to the phase variation 
between eclipses using the model given in equation (f). These 
results are for {N x , N v , N xy ) = (5, 10, 5) with positions measured 
using the CNTRD method and an aperture radius of 2.25 pixels. N 
is the number of points included in the fit and BIC is the Bayesian 
information criterion as defined in the text. A is the amplitude 
of the thermal phase effect and </>max is the phase of maximum 
brightness relative to phase 0.5. Random and systematic errors 
are given for each quantity in that order. 



Parameter 


Channel 1 


Channel 2 


ai [%] 


0.148 ±0.005 ±0.01 


0.183 ±0.004 ±0.01 


a 2 [%] 


0.003 ±0.005 ±0.01 


0.023 ±0.005 ±0.01 


bi [%] 


0.001 ±0.003 ±0.01 


-0.006 ±0.003 ±0.01 


A [%] 


0.296 ±0.009 ±0.02 


0.366 ±0.007 ±0.02 


^raax 


0.001 ±0.003 ±0.01 


-0.010 ±0.006 ±0.02 


JV 


133124 


179851 


x 2 


142398.6 


216717.3 


BIC 


142928.2 


217259.0 


RMS [%] 


0.539 


0.717 



aperture radius and position measurement method for our 
various lightcurves. Also plotted in Fig. [5] are the ampli- 
tude of the phase variation and the offset from phase 0.5 
to the time of maximum brightness in phase units. There 
is some dependence on aperture radius for these results, 
e.g., the value of ax for both channels show a trend to- 
wards smaller values with increasing radius. We also see 
that there is worse agreement between the results for dif- 
ferent position measurement methods for smaller apertures 
as a result of the increased sensitivity of the pixelation noise 
to small differences in the assumed position. For all of the 
coefficients in both channels we see that the results vary by 
about ±0.01 per cent as a function of aperture radius. We 
therefore assume that systematic noise limits the accuracy 
of these results to ±0.01 percent. 

Despite the limit of ±0.01 per cent in the accuracy of 
these results, we are able to draw some definite conclusions 
about the thermal phase effect in WASP-18. Firstly, the am- 
plitude of the thermal phase effect is very similar to the 
depth of the secondary eclipse. This can be seen in Fig. 2] 
and by comparing the values for the amplitude in Table [2] 
to the eclipse depths given in Table [3] Secondly, the off- 
set between phase 0.5 and the time of maximum brightness 
due to the thermal phase effect is consistent with to within 
about 0.01 phase units for the channel 1 data and 0.02 phase 
units for the channel 2 data. Thirdly, the parameter a 2 is 
also consistent with the value so the shape of the thermal 
phase variation is sinusoidal to within the limits sets by the 
systematic noise. 



3.9 Eclipse model 

We tried several different methods to model the entire 
lightcurve for each channel including both eclipses, the phase 
variation and PDSV, but were not able to find any method 
that gave reliable results. We suspect that there is some fac- 
tor other than position on the detector that introduces sys- 
tematic noise at the level of ~0.01 percent with a timescale 
of ~day. This can be seen in Fig. [4] where there are clear 



systematic errors remaining in the lightcurve at this level. 
These systematic errors are not removed by increasing the 
complexity of the model used for the PDSV. It will be diffi- 
cult to identify this additional factor given that little infor- 
mation about the shape of the point spread function can be 
measured from the undersampled IRAC images. 

There are many published studies that have used IRAC 
photometry obtained over ~ 5 hours of observation to suc- 
cessfully model hot Jupiter eclipses, so we decided to only 
model the data within 0.1 phase units of the primary and 
secondary eclipses. We fit these data simultaneously using 
a single model to account for the true flux variations of the 
system. We then account for systematic errors in the mea- 
sured flux independently for the data around each eclipse. 

We have used the NDE lightcurve mo del 
l|Nelson fc Davisll 19721 ; Fopper fc Etzelll98ll ; lEtze]||l98lh to 
model the primary and secondary eclipses in our lightcurves. 
This model uses biaxial ellips oids to approxim ate to pro- 
jected area of the star/planet. iGimenea (|2006T ) has shown 
that this model used with an integration ring size of one 
degree (as we have done) can be used to model planetary 
transits with a precision of ~ 4 x 10 -5 , which is sufficient 
for our purposes. From inspection of our model lightcurves 
we find that ~ 4 of the model data points during primary 
eclipse are affected by numerical noise at this level and that 
there is no numerical noise during secondary eclipse. We 
created a double-precision version of the NDE model that 
has negligible numerical noise, but that runs appreciably 
slower than the original single-precision code. We used our 
double-precision version to verify that the numerical noise 
in the single-precision version has a negligible effect on our 
results, so all the results presented here are based on the 
single-precision version. 

We use the NDE model to calculate i s (cj>) and £ p (4>), the 
contribution of the star and planet, respectively, to the total 
apparent flux at any given phase, <j>, including the effects of 
tidal distortion and eclipses. Note that l s and £ p include 
the effects of the eclipses and transits and the ellipsoidal 
variation of both star and planet. To model the variation in 
magnitude due to the phase effect of the planet we use the 
harmonic series 

h(<j>) = ai cos(0) + bi sm(cj>) + a 2 cos(2^) — /i max , 

where the values of ai, 6i and a 2 are taken from the least- 
squares fit to the data between transit and eclipse for the 
same aperture size and position measurement method and 
/imax is chosen such that the maximum value (corresponding 
to minimum flux) of h(<j>) is 0. 

The apparent magnitude of the system is then given by 

m, = m - 2.5 log [4(^0 + tp(<f>i)] + ip((t>i)h((()i)/ip,max, 

where f p , ma x is a normalization factor. Our calculations are 
done using magnitudes but we present the results in flux 
units or as percentages. 

In addition, we model the PDSV independently for the 
data around primary and secondary eclipse using Legendre 
polynomial functions of the x and y position plus an optional 
linear function of time. For each set of lightcurve model 
parameters we calculate the optimum values of the PDSV 
model parameters using singular value decomposition to fit 
the residuals from the lightcurve model. 

The parameters of the NDE lightcurve model of rele- 
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Figure 4. The flux of WASP-18 measured in IRAC channel 1 (left panel) and channel 2 (right panel) in an aperture of radius 2.25 pixels 
after correction for position dependent sensitivity variations (PDSV) for the parameter sets (7V cos , -/V a ; n , N x , N y , N xy ) = (2,1,5,10,5) 
and positions measured with the CNTRD method. Data are plotted averaged in 200 s bins for clarity and the best-fit sinusoidal model is 
also shown. The mean value in secondary eclipse is indicated with a dotted line. Note that data in eclipse (small points) are not included 
in the fit. The PDSV model is shown as a function of time in the middle panels and as a function of position as a grey-scale plot in the 
lower panels. The grey-scale is linear between ±4 per cent for channel 1 and 1 per cent for channel 2 with positive values being white. 



vance to our study are: J, the surface brightness of the planet 
in units of the central surface brightness of the star exclud- 
ing the thermal phase contribution; r± — i? a tar/a, the radius 
of the star in units of the semi-major axis; T2 = 7? p ianot/a, 
the radius of the planet in units of the semi-major axis; i, 
the inclination, u*, the linear limb-darkening coefficient for 
the star; ecos(o;) and esin(oj), where e is the orbital ec- 
centricity and lj is the longitude of periastron. We fix the 
mass ratio of the system at the value m p i ane t /wi s tar = 0.01. 
We did not use this combination of parameters directly as 
free parameters in our least-squares fitting because there are 
significant correlations between them. Instead, we introduce 



the following parameters which are more directly related to 
the observed features of the lightcurve: Am tr , Am oc , W, 
S. The parameters of the NDE lightcurve model are then 
calculated as follows: 

r-i = fcn; 

, _ /(l-fc) 2 -S 2 (l + fc) 2 

b - V i - ^ ; 
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Figure 5. Coefficients of the sinusoidal model for the phase variation in WASP-18 as a function of aperture radius. The parameter 
set (N x , N y , N X y) = (5,10,5) was used for the correction for position dependent sensitivity. Plotting symbols are as follows: CNTRD — 
squares; GCNTRD - diamonds; GAUSS2D - filled circles. The solid line in the upper panels shows the semi-amplitude of the phase variation 
for the CNTRD results. The solid line in the lower panel is the the phase offset from phase 0.5 for time of maximum brightness for the 
thermal phase effect derived from the CNTRD results. Points have been offset horizontally by ±0.05 for clarity. 



i = cos 1 (6ri) ; 
j _ 1 - «*/3 

1.2 ( 2.5 _ i\ 

I ln(10)Am oc X J 

These parameters are adapted from 

ISeager fc Mallen-Ornelasl (|2003h so that, for a circular 
orbit, Amtr is the depth of the primary eclipse in mag- 
nitudes, Am oc is magnitude difference between the flux 
duration occultation and the minimum of the thermal 
phase curve; W is the width of the transit in phase units 
and 5* is the duration of the ingress phase of the transit in 
units of W. The intermediate variables used here are k, the 
radius ratio and b, the impact parameter. We also include 
a correction to the time of mid-transit compared to the 
ephemeris of Hellier et al., A2~b. 

We are careful here to define what we mean by the depth 
of the secondary eclipse because the variation in flux due 
to the thermal phase effect on the timescale of the eclipse 
is comparable to the precision with which we can measure 
the depth from our photometry (~ 0.01 per cent) and the 
maximum of the thermal phase effect may not occur at 
mid-eclipse. For ease of calculation, interpretation and com- 



parison with other measurements, we simply measure the 
mean flux during occultation (excluding ingress and egress 
phases), /i n , and the mean flux either side of the eclipse in a 
region ± 0.1 phase units around the time of mid-eclipse, / ou t 
and define the secondary eclipse depth to be D = £aat=zi£ . 

/out 

The fluxes are measured from the lightcurve corrected for 
PDSV. 

There are some second-order effects not accounted for 
by our model. We do not account for the brightness distri- 
bution on the day-side of the planet, but this will have a 
negligible effect on our results given that the thermal phase 
effect is not strongly peaked and is symmetrical about phase 
0.5, so this distribution will be approximately uniform and 
symmetrical. Doppler boosting is negligible compared to our 
signal-to-noise (< 0.001 per cent). We make a small correc- 
tion to the results for the effects of image persistence in the 
channel 2 data by assuming a dilution of the eclipses for 1±1 
percent. We have applied a correction to the apparent times 
of secondary eclipse for the light travel time across the or- 
bit of 2a/ c ~ 20 s so that the times and phases quoted here 
are the true time of mid-occultation relative to the apparent 
time of mid-transit. We also apply a correction to the values 
of ecos(oj) quoted below for this light travel time. 
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Figure 6. Upper panel: Raw photometry for an aperture radius of 3 pixels (filled symbols - channel 1, open symbols - channel 2) 
together with the correction for PDSV based on positions measured using the CNTRD method (lines). Lower panel: Photometry corrected 
for PDSV (points) and models fit by least-squares (lines). The channel 2 data have been vertically offset by 0.03 in the upper panel and 
0.01 in the lower panel. In both panels the data and models are plotted in 60s bins. 



We use the simplex algorithm of lNelder fc Meadl (|l965l ) 
to optimise the least-squares fit of our model to the 
lightcurves. The simplex algorithm is a simple way to op- 
timize a least-squares solution given an initial set of param- 
eters, but it is not guaranteed to find the global minimum 
value of x 2 in the parameter space. In this case we are able 
to estimate accurate initial values for the most important 
parameters and so any solution will not be very far from the 
global minimum. However, we do find that numerical noise 
prevents us from using this algorithm by itself to find the 
optimum solution. We work around this problem by test- 
ing many initial starting values. We found that the solution 
with the lowest value of \ 2 sometimes has parameters that 
are slightly biased when compared to other solutions with 
similar values of \ 2 as a result of the numerical noise. We 
avoid this problem by taking the median value of each pa- 
rameter for all solutions with \ 2 within 5 of the minimum 
as our best estimate of the parameter. 

The results for the depths of the eclipses and the RMS 
of the residuals for each data set are shown in Fig. [7] The 
depth of the primary eclipse (transit) can vary slightly with 
wavelength because the apparent radius of the planet will 
be large r at wavelengths where t he atmosphere has a large 
opacity l|Seager fc Sasselovl f2000). The size of this effect is 
approximately 2//7? p i anct /J?s tar where H is the atmospheric 
scale height. In practice, for WASP-18b the size of this ef- 



fect is negligible (~ 0.001 percent) because the large surface 
gravity of this massive planet makes the scale height of the 
atmosphere (~ 40 km) much smaller than the size of the star. 
It can be seen from Fig.[7|that the solutions with the lowest 
RMS occur for an aperture radius 2 pixels, but the transit 
depths for channel 1 and 2 disagree by about 0.01 per cent 
for these data sets. 

The transit depths measured in channels 1 and 2 are 
consistent with each other for an aperture radius of 3 pixels 
and lie near the centre of the range of values obtained. The 
best fit to the lightcurves for an aperture radius of 3 pixels 
using the positions from the GCNTRD method are shown in 
Fig. [6] and the parameters for the model used are given in 
Table [3] It can be seen that there is some residual corre- 
lated noise in the lightcurves after removal of the model 
for the PDSV. We quantified this residual correlated noise 
by calculating the RMS of the resid uals after binning for a 
range of bin sizes (|Pont et alj|200fj ). The results are shown 
in Fig. [8] and compared to the expectation for pure photon 
noise. At the timescale of the eclipse it can be seen that the 
channel 2 data are only weakly affected by correlated noise 
(< 0.005 per cent) but the channel 1 data are affected corre- 
lated noise at a level of 0.005-0.01 percent, particularly for 
the data covering occultation. 

Given that there is significant correlated noise in our 
data, we decided to calculate the random error on our model 
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Figure 7. Depths of eclipses for transit and occultation mea- 
sured by fitting the eclipses. Channel 1 data are shown with filled 
symbols, channel 2 data with open symbols. Different symbols 
denote different position measurement methods. The RMS of the 
residuals of the fits are also shown using the same symbols. 



param eters using the "prayer-bead" method (jPont et al.l 
2006). This uses a circular permutation of the residuals by 
a random number of steps to create mock data sets. We 
applied the circular permutation to the residuals of the pri- 
mary and secondary eclipses independently and then used 
the simplex algorithm to fit models to 1024 mock data sets. 
The standard deviation of parameters from the fits is used 
to calculate the random errors for the model parameters 
given in Table [3] based on the analysis of the lightcurves for 
an aperture radius of 3 pixels using the positions from the 
GCNTRD method. The random errors quoted include the ef- 
fect of the uncertainty in correcting for image persistence 
in the data. We use the range of values from the different 
apertures and position measurement methods to estimate 
the systematic errors on each parameter. The distribution 
of the parameters for the mock data sets is shown for some 
parameters of interest in Fig. [9] As can be seen from this fig- 
ure, the eclipse depths derived from the mock data sets can 
be biased by up to ~ 0.005 per cent from the actual value. 
This is an consequence of the correlated noise in the resid- 
uals. We also used the results from these mock data sets to 
calculate the Pearson correlation coefficient, r, for all pairs 
of free parameters used in the least-squares fit. There is a 
weak anti-correlation between W and S (r « —0.5), and a 
weak correlation between W and esin(cj) (r ~ 0.5), but the 
other free parameters are uncorrelated, as expected. 



4 OPTICAL VARIABILITY 

The interpretation of our data would be considerably com- 
plicated by any intrinsic variability of the star WASP-18. 

We have analysed the WASP lightcurves of WASP-18 
to determine whether they show periodic modulation due 
to the combination of magnetic activity and the rotation 
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Figure 8. RMS of the residuals after binning as a function of bin 
size (solid line) compared to the predicted photon noise (dotted 
line). The vertical, dashed lined in each panel shows the duration 
of eclipse and the duration of ingress/egress. 
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Figure 9. Parameter correlation plots from our residual permuta- 
tion error analysis for channel 1 (blue crosses) and channel 2 (red 
diamonds). Our adopted values for each channel are indicated 
using dotted and dashed lines for channel 1 and 2, respectively. 



of the star. The o bserved value of V ro t sin J = llkms 1 
ijHellier et alJ l2009) together with the stellar radius imply a 
rotation period of about 6 d for WASP-18. We used t he sine- 
wave fitting method described in lMaxted et alj |201ll ') to cal- 
culate periodograms over 4096 uniformly spaced frequencies 
from to 1.5 cycles/day. The false alarm probability (FAP) 
for the strongest peak in these periodograms was calculated 
using a boot- s trap Monte Carlo method also described in 
iMaxted et ail (^Olll ). Variability due to star spots is not ex- 
pected to be coherent on long timescales as a consequence 
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Table 3. Results of least squares fits to the primary and secondary eclipses. J' = J /(I — m*/3) is the ratio of the integrated surface 
brightness of the star and the day-side of the planet. Other symbols are defined in the text. Random and systematic errors are given for 
each parameter in that order. Parameters that can be derived from the analysis of the optical TRAPPIST lightcurves are also given in 
the final column. 



Parameter 


Channel 1 


Channel 2 


TRAPPIST 


Am tr [%] 


0.969 ±0.013 ±0.007 


0.979 ± 0.013 ± 0.009 


0.965 ± 0.056 


Am oc [%] 


0.015 ± 0.014 ± 0.009 


0.028 ±0.006 ±0.017 




D [/o] 


n 3fM 4- n m t 4- n nno 
U.OU4 it U.U1 1 ± U.UUy 


n qvo 4- n nna 4_ n m q 
U.oiy it U.UUo it U.Ulo 




1 I 


u.uyoo It U.UUUD It U.UUUo 


n nO/io 4- n nnnQ 4- n nnno 
u.uy4z it u.uuuo it u.uuuz 


U.Uy40 it U.UUll 


c 
O 


U. (92 ± 0.U09 ± U.OUO/ 


0.8U2 ± U.OUo ± 0.U0/ 




u* 


0.06 ±0.03 ±0.06 


0.07 ±0.03 ±0.04 




e cos (a;) 


0.0002 ± 0.0004 ± 0.0003 


0.0001 ± 0.0002 ± 0.0003 




e sin(w) 


-0.003 ± 0.006 ± 0.004 


-0.001 ±0.003 ±0.002 




AT [s] 


-109±8±0 


-108±8±0 




ri 


0.287 ± 0.006 ± 0.004 


0.282 ± 0.005 ± 0.004 


0.291 ±0.017 


V2 


0.0281 ± 0.0006 ± 0.0003 


0.0278 ± 0.0005 ± 0.0005 


0.0286 ± 0.0026 


k 


0.0982 ± 0.0007 ± 0.0004 


0.0987 ± 0.0006 ± 0.0004 


0.0983 ± 0.0030 


b 


0.39 ±0.06 ±0.04 


0.32 ±0.06 ±0.04 


0.41 ± 0.15 


i 


83.6 ±1.0 ±0.7 


84.8 ±1.1 ±0.8 


83 ±3 


J' 


0.33 ±0.02 ±0.02 


0.40 ±0.01 ±0.02 




Phase of mid-occultation 


0.5003 ± 0.0006 ± 0.0004 


0.5002 ± 0.0003 ± 0.0005 




^transit 


33996 


35261 


3649 


-^occult at ion 


36847 


34355 




x 2 


74213.9 


90739.5 


2768.5 


RMS [%] 


0.55 


0.75 


0.40 



of the finite lifetime of star-spots and differential rotation 
in the photosphere and so we analysed the data from each 
observing season independently. We removed the transit sig- 
nal from the data prior to calculating the periodograms by 
subtracting a simple transit model fro m the lightcu r ve. In 
addition to the 2 seasons of data from iHellier et ail (^OOOt ) 
we also analyse 6041 observations obtained during the 2012 
Jun - Dec observing season. This date range covers the time 
of our Spitzer channel 2 observations. 

We did not find any significant periodic signals (FAP< 
0.05) for WASP-18 apart from frequencies near 1 cycle/day 
due to instrumental effects. We examined the distribution 
of amplitudes for the most significant frequency in each 
Monte Carlo trial and used these results to estimate a 
95 per cent upper confidence limit of 0.1 per cent for the am- 
plitude of any period ic signal in these WASP lightcurves. 
iBeaulieu etall l|201uft have shown that the amplitude of the 
modulation at IRAC wavelengths due to star spots in solar- 
type stars is an order-of-magnitude smaller than at optical 
wavelengths. We conclude that any intrinsic variability of 
WASP-18 due to star spots has a negligible impact on our 
analysis. 



5 ECLIPSE EPHEMERIS 

The analysis above provided two new, precise measurements 
of the time of mid-occultation and mid-transit. In addition 
we have 5 new times of mid-transit from our TRAPPIST 
observations. A global analysis of the five lightcurves was 
performed with the MCMC software described by Gillon et 
al. (2012). In addition to the baseline model and to the tran- 
sit ephemeris and shape parameters, the timings of the tran- 
sits were included as free parameters, the transit ephemeris 



being constr ained by norma l prior s based on the ephemeris 
presented bv lNvmever et all (|201ll ). Th e details of th i s ana l- 
ysis are similar to the ones described in iGillon et al.l (|2012h . 
The parameters derived from the least-squares fit to the 5 
lightcurves are shown in Table [3] It can be seen that there is 
very good agreement between the parameters of the system 
dervied from the optical and infrared lightcurves. 

We have also measured a new time of mid-transit by 
using a least-squares fit of the NDE lightcurve model to the 
2010 season of WASP data. The time of mid-transit quoted 
is close to the mid-point of dates of observation for these 
data and the standard error on the time of minimum is cal- 
culated using the prayer-bead method. All these times of 
mid-eclipse are given in Table [4] together with other pub- 
lished times of mid-eclipse. 

We used a least squares fit with a single value of the 
period and the times of mid-transit and mid-occultation as 
free parameters to determine a following linear ephemeris. 

TDB(mid - transit) = 2455265.5525(1) ± 0.9414523(3) • E 
TDB(mid - occult.) = 2455266.0234(3) + 0.9414523(3) ■ E 



The x 2 value for this fit was 21.8 with 11 degrees-of- 
freedom so the standard errors q uoted in the final digits here 
have been scaled by y2L8/ll. We also tried a quadratic 
ephemeris fit to the same data but found that this did not 
significantly improve the fit. The residuals from this O — C 
diagram for this linear ephemeris is shown in Fig. 1101 It 
can be seen that the times of eclipse for WASP-18 have not 
varied by more than about 100 s over 3 years. 
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Table 4. Apparent Barycentric Dynamical Time (TDB) of mid- 
transits (tr) and mid-occultation (oc) for WASP-18. The cycle 
number is calculated from our updated linear cphemeris and O—C 
is the residual from this ephemeris. Times of mid-occultation have 
been corrected for light- travel time across the orbit. 



BJDtdb 


- 2450000 


Type 


Cycle 


O-C 


Source 


4664.9061 


±0.0002 


tr 


-576 


0.00013 


1 


4820.7168 


±0.0007 


oc 


-411 


0.00019 


2 


4824.4815 


±0.0006 


oc 


-407 


-0.00097 


2 


5220.8337 


±0.0006 


oc 


14 


-0.00006 


a 


5221.3042 


±0.0001 


tr 


15 


0.00017 


a 


5392.6474 


±0.0004 


tr 


197 


0.00014 


3 


5419.0083 


±0.0012 


tr 


225 


0.00015 


b 


5431.7191 


±0.0003 


oc 


238 


-0.00114 


a 


5432.1897 


±0.0001 


tr 


239 


-0.00087 


a 


5470.7885 


±0.0004 


tr 


280 


-0.00057 


c 


5473.6144 


±0.0009 


tr 


283 


0.00098 


c 


5554.5786 


±0.0005 


tr 


369 


0.00026 


c 


5570.5842 


±0.0006 


tr 


386 


0.00118 


c 


5876.5559 


±0.0013 


tr 


711 


0.00097 


c 
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Figure 10. Residuals from our best-fit linear ephemeris for ob- 
served times of mid-transit (circles) and mid-occultation (squares) 
for WASP-18. The difference between our ephemeris and the 
ephemeris of Hellier et al. is also shown (dashed line). 



6 WASP-18 STELLAR PARAMETERS 

A total of 21 individual HARPS spectra of WASP-18 were 
co-added to produce a single spectrum with a typical S/N of 
around 200: 1. The analy s is wa s performed using the meth- 
ods given in lDovle et al.l (|2012l ). The H a and Hp lines were 
used to give an initial estimate of the effective temperature 
(T e ff). The surface gravity (logg ) was determined fr om the 
Ca I lines at 6162A and 6439A (|Bruntt et al.ll2010r i, along 
with the Na I D lines. Additional T e ff and logt/ diagnostics 
were performed using the Fe lines. An ionisation balance 
between Fe I and Fe n was required, along with a null de- 
pendence of the abu ndance on either e quivalent width or 
excitation potential ijBruntt et al.l[2008l ). This null depen- 
dence was also required to determine the micro-turbulence 
(£t). The parameters obtained from the analysis are listed 



Table 5. Stellar parameters of WASP-18 from our spectroscopic 
analysis. 



Parameter 




Value 




[K] 


6400 ± 75 


logg 




4.29 ± 0.10 


£t 


[kms -1 ] 


1.20 ± 0.08 


v sin i 


[kms -1 ] 


12.1 ± 0.5 


[Fe/H] a 




0.10 ± 0.08 


Mass b 


[Mq] 


1.26 ± 0.09 


Radius' 


[Rq] 


1.25 ± 0.15 


Sp. Type c 


F6 


Distance 


[PC] 


130 ± 20 



a [Fe/H ] is relative to the solar value obtaine d by lAsp lund e t al.l 
(l2009h . b Mass and radius estimated using the lTorres et al.l ll20ldl ) 
ca libration. c Spectral Type estimated from T e g using the table 
in lGravl fcOOSl) . 



in Table [S] The value of [Fe/H] was determined from equiv- 
alent width measurements of several unblended lines, and 
additional least squares fitting of lines was performed when 
required. The quoted error estimates include that given by 
the uncertainties in T c s, logg, and £ t , as well as the scatter 
due to measurement and atomic data uncertainties. 

The projected stellar rotation velocity (v sin i) was de- 
termined by fitting the profiles of several unblended Fe I 
lines. A value for macroturbulence (iWc) of 4.6 ± 0.3 kms" 1 
was a ssumed, based on the calibration by iBruntt et alj 
|2010T) . An instrumental FWHM of 0.07 ± 0.01 A was deter- 
mined from the telluric lines around 6300A. A best fitting 
value of vsini = 10.9 ± 0.7 kms -1 was obtained. 

The rotation rate (P — 5.8 ± 0.8 d) implied by the 
vsini g ives a gyrochr onological age of ~ 1.1 Jg -7 Gyr us- 



ing the H arnc sf (|2007D relation. The value of T e ff derived 
from our spectroscopi c analysis agrees we ll with the value 
6455±70 K derived bv lMaxted et al.1 l|201ll ) from optical and 
near-infrared photometry using the infrared flux method. 
The distance derived here assuming that WASP-18 is a 
main-sequence star and quoted in Table [5] is consistent with 
the value 100 ± 10 pc derived from the Hipparcos parallax 
(van Leeuwenl [20071 ). 



7 PHYSICAL PARAMETERS 

Our new photometric and spectroscopic results allow for 
an improved determination of the physical properties of 
the WASP-18 syste m. We performed this analysis follow- 
ing the method of ISouthworthl ((2009) , which requires as 
its input parameters measured from the lightcurves and 
spectra, plus tabulated predictions of theoretical models. 
From the lightcurves we adopted ri = 0.284 ± 0.005, r-z = 
0.0280 ±0.0005 and i = 84° ± 1°. The stellar T cB and [Fe/H] 
were taken from the spectroscopic determination in the pre- 
vious section, and the star's v elocity amplitude w as taken 
to be K x = 1816.7 ± 1.9 ms -1 (|Triaud et al.ll201fj| ). 

An initial value of the velocity amplitude of the planet, 
K-2 , was used to calculate the physical p roperties of the sys - 
tem with the physical constants listed bv lSouthworthl l|201ll ). 
The mass and [Fe/H] value of the star were then used to ob- 
tain the expected T c ff and radius, by interpolation within 
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Table 6. Derived physical properties of the WASP-18 system. 
Parameter values are shown with random and, where appropriate, 
systematic errors, respectively. 



Parameter 


Value 


M A (Mq) 


1.295 ± 0.052 ±0.027 


Ra (Rq) 


1.255 ± 0.027 ±0.009 


log 9a (cgs) 


4.353 ± 0.017 ±0.003 


PA (p©) 


0.655 ±0.035 


M b (M Jup ) 


10.52 ± 0.28 ±0.15 


Rh (-Rjup) 
9h (ms~ 2 ) 


1.204 ± 0.027 ±0.008 


179.9 ± 6.4 


Pb (P.Jup) 


5.64 ± 0.31 ±0.04 


r cq (K) 


2411 ± 35 


a (AU) 


0.02055 ± 0.00028 ± 0.00014 


Age (Gyr) 


n A +0.8 +0.5 
-0.9 -0.3 



one set of tabulated predictions from stellar theory. Ka was 
refined in order to find the best agreement between the ob- 
served and expected T e g, and the measured n and expected 
— . This was performed for ages ranging from the zero-age 
main sequence to when the star was significantly evolved 
(log g < 3.5), in steps of 0.01 Gyr. The overall best fit was 
found, yielding estimates of the system parameters and also 
the stellar age. 

This procedure was performed separately u sing five dif- 
feren t sets of stellar theoretical models (see ISouthworthl 
|2010| . for details) plus a calibration of stellar prop- 
erties based on well-st udied eclipsing binary star sys- 
tems llEnoch et alllioioh . with calibration coefficients from 
ISouthworthl (|201ll ). The results are given in Tabled where 
we quote the mean value for each parameter, the random er- 
ror and an estimate of the systematic error from the range 
of values derived from the different stellar models, where ap- 
propriate. It can be seen from Table [B] that the results from 
different models are consistent to within the random errors 
on each parameter. 

In comparison to previous work, we have derived more 
precise radii, surface gravities and densities, for both com- 
ponents. We constrain the age of the star to be less than 
1.7 Gyr, consistent with the gyrochronological age derived 
above. 



8 POSSIBILITY OF CONTAMINATION BY A 
COMPANION STAR 

We have estimated the probability that our Spitzer photom- 
etry of WASP-18 is contaminated by the "third-light" from a 
companion star. It is not possible to detect a modest amount 
of third-light contamination directly from the lightcurve it- 
self because its only effect is to reduce the depths of the 
eclipses. It would be possible to find a good fit to a lightcurve 
affected by third-light contamination, but the parameters of 
the lightcurve model would be biased, e.g., k would be too 
small. 

Our calculation is based on the upper limit from our 
AO observations of 4.0 magnitudes for the brightness of any 
companion between 0.2-2arcseconds from WASP-18 and 
the upper limit of 43 ms -1 y _1 over a baseline of 500 days 
to the variation in the mean radial velocity of WASP-18 



from iTriaud et af] (|20ld ). We assume that the probability 
distribution for the mass, eccentricity and period of the hy- 
pothesised companion is the same as the distributions fo r 
companions to solar- type stars from iRaghavan et al.ll|2010h . 
We approximated the distribution of companion masses us- 
ing a uniform distribution from 0.2 Mq to 0.8 Mq and used 
a uniform eccentricity distribution from to 1. We then 
created a set of 65536 simulated binary stars with randomly 
selected periods, masses and eccentricities according to these 
distributions and randomly orientated orbits. We found that 
of these simulated binary stars, approximately 55 per cent 
would have been resolved by our AO imaging at the dis- 
tance of WASP-18; 20 percent would have orbital periods 
less than 500 days and a semi-amplitude of 43ms -1 or more 
and 25 per cent would show a change in radial velocity of 
43ms -1 or more over 500 days. This leaves only 5 per cent 
of hypothesised binaries as not detectable given our AO 
imaging and the published radial velocity data. The proba- 
bility that WASP-18 has a stellar companion is further re- 
duced because the overall binary fraction ob served for planet 
hosti ng stars is approximately 25 per cent (|Raghavan et al.l 
120101 ). An M-dwarf at the same distance as WASP-18 just 
below out detection limit of 4.0 magnitudes in the K-band 
at 0.2 arcsec would contribute no more than 5 per cent of the 
light at 4.5 /im. The more stringent limit of 6.0 magnitudes 
in the K-band that applies for separations of 0.5 -2.0 arcsec 
corresponds to an M-dwarf that contributes no more than 
1 per cent at 4.5 pm. Of the simulated binary stars approxi- 
mately 45 per cent would be detected at this resolution. 

In conclusion, our AO imaging and the published ra- 
dial velocity data show that it is unlikely that WASP-18 
has a stellar companion that significantly contaminates our 
Spitzer photometry. 



9 DISCUSSION 

The values of D in Table [3] are in very good agreement 
the values 0.31 ± 0.02 at 3.6 /zm an d 0.38±0. 02 percent a t 
4.5 /dm measured independently by iNvmever et all (|201lh . 
Their analysis of the secondary eclipse depths in 4 IRAC 
passbands suggests that the day-side atmosphere of WASP- 
18 is likely to feature a temperature inversion. For zero 
albedo and zero redistribution of heat to the night side of the 
planet the integrated brightness temperature for the day- 
side is T E=0 = (§) 1/4 To = 3110 ± 35 K (|Cowan fc Agoll 
2011). For black-body emission this implies eclipse depths 
of 0.329 ± 0.005 at 3.6 ^m and 0.379 ± 0.011 per cent at 
4.5 /jm, both in good agreement with the observed values. 
Zero-redistribution of heat within the atmosphere is also 
consistent with our observation that the peak of the ther- 
mal phase curve occurs close to the time of mid-occultation. 
Little can be said about the chemical composition of the 
day-side atmosphere at this stage because no strong molec- 
ular absorption or emission features have been detected from 
these secondary eclipse depth measurements. 

The amplitude of the thermal phase curve we have 
measured and the lack of a significant offset between the 
maximum in this curve and the time of mid-eclipse are 
both consistent with the conclusion based on the secondary 
eclipse depths that the albedo and recirculation efficiency 
for WASP-18 are both very low. This is consistent with 
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the hypothesis that very hot Jupiters have weak recircula- 
tion b ased mainly on second ary eclipse depth measurements 
alone (|Cowan fe Agol|[201ll ). The good agreement between 
the recirculation efficiency inferred from the eclipse depths 
and measured from the thermal phase curve for WASP-18 
strengthens this conclusion. 

The stellar limb darkening at infrared wavelengths 
is lower than at optical wavelengths and so the tran- 
sit produces a more "box-shaped" eclipse. This, combined 
with the precise photometry that is possible with Spitzer 
IRAC data, results in more precise estimates for param- 
eters SUch as Rstar/d, iiplanct/a and k = i?planct/fl!star. 

Our res ults for these paramete rs agree well with the re- 
sults of ISouthworth et al.l (120091 ). The agreement with the 
results of iTriaud et al.l ([2010) is less good mainly because 
they find a larger stellar radius than our study (R s tm/a = 
0.313 ± 0.010). The values of ecos(oj) an d esin(u) der ived 
from our analysis agree well with those of lTriaud et al.l . but 
the value of esin(oj) = 0.0085 ± 0.0009 they derive from 
their high quality radial velocity data is much more precise 
than ours a nd points to a sma ll but significantly non-zero 
eccentricity. lArras et all l|2012h have argued that the small 
value of this apparent eccentricity combined with longitude 
of periastron very close to uj = 90° is exactly the signal ex- 
pected due to surface flows induced by tides on the planet. 
Their conclusion that the orbital eccentricity of WASP-18 b 
is less than 0.009 is consistent with the results of our analy- 
sis, although we are not able to confirm whether e <C 0.009 
as they suggest. 

The measurement of the thermal phase effect for hot 
Jupiters using a continuous set of observations over an or- 
bital cycle with Warm Spitzer is not a well-established tech- 
nique, so it is useful to compare our experience of observing 
WASP-18 with the results usi ng a similar observi ng strat- 
egy for WASP - 12 obtained by I Cowan et all (|2012h and for 
HD 189733 bv lKnutson et all (|2012T l. We find that system- 
atic errors of unknown origin limit the accuracy with which 
we can measure the amplitude of signals with time-scale 
comparable to the orbi tal period to about ±0.01 percent. 
The main difficulty that lCowan et all report in their WASP- 
12 analysis is a signal on twice the orbital frequency in 
the 4.5 /im data that they tentatively attribute to the el- 
lipsoidal modulation of WASP-12b. However, as they make 
clear, this signal is not seen in their 3.6 /jm data and 
may be due to "uncorrected systematic noise". The ampli- 
tude of this "cos(2</>)" signal in their 4.5 /jm data is about 
±0.1 percent, ten times larger than the level of systematic 
noise we find on these timescales. We do not see any sig- 
nal for the cos(20) harm onic in our 4.5 /jm data greater 
than about 0.02 per cent. [ Knutson et alj use the wavelet- 
based method of ICarter fc Winnl |2009) to account for sys- 
tematic noise in their full-orbit lightcurves of HD 189733 by 
assuming that this noise has a power spectral density vary- 
ing as 1/ frequency. They find that the systematic noise con- 
tributes 0.0162 percent of the total scatter in their channel 
1 data - comparable to the level seen in our data - but 
only 0.0017 per cent in channel 2 - much less than we see 
in our data. A systematic application of the wavelet-based 
method to archival Spitzer data may be a useful way to 
better understand the systematic noise sources in this in- 
strument and perhaps identify observing strategies that can 
reduce systematic noise levels. 



10 CONCLUSIONS 

The amplitude, shape and phase of the thermal phase ef- 
fect we have measured from our Warm Spitzer lightcurves 
of WASP-18 are consistent with a sinusoidal variation with 
the same amplitude as and symmetric about the secondary 
eclipse, to within an accuracy ~ 0.01 per cent set by some 
unknown source of systematic error. One contribution to 
this systematic error is likely to be the image persistence we 
observe from the offset images obtained immediately after 
our observations of WASP-18. This leads to the conclusion 
that WASP-18 b has a low albedo and that heat transport 
to the night-side of the planet is ineffic i ent. T his is the same 
conclusion reached by iNvmever et al.l (|201lf ) based on the 
eclipse depths at 3.6 /jm, 4.5 /im, 5.8 /im and 8.0 fim. The 
eclipse depths we measure at 3.6 /im and 4.5 fim are consis- 
tent with the previous measurements. 
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